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Supersymmetry predicts the existence of supersymmetric partners for each of the par- 
ticles of the standard model. The search for these sparticles is a principal motivation of 
the forthcoming Run II of the Fermilab Tevatron collider and of the CERN Large Hadron 
Collider (LHC) program. A potentially important, but heretofore largely overlooked, dis- 
covery channel is the associated production of a spin-1/2 gaugino (x) with a spin- 1/2 gluino 
(g) or with a spin-0 squark (q). Color- neutral gauginos couple with electroweak strength, 
whereas the colored squarks and gluinos couple strongly. Associated production is therefore 
a semi-weak process in that it involves one somewhat smaller coupling constant than the pair 
production of colored sparticles. However, in popular models of SUSY breaking ]l]||, the 
mass spectrum favors much lighter masses for the low-lying neutralinos and charginos than 
for the squarks and gluinos. This mass hierarchy means that the phase space for production 
of neutralinos and charginos, the corresponding partonic luminosities, and the production 
cross sections will be greater than those for gluinos and squarks. These advantages are 
potentially decisive at a collider with limited energy, such as the Tevatron. Furthermore, 
associated production has a clean experimental signature. For example, the lowest lying 
neutralino is the (stable) lightest supersymmetric particle (LSP) in supergravity (SUGRA) 
models [TJ, manifest as missing energy in the events, and it is the second lightest in gauge- 
mediated models f2|. In models with a very light gluino ||, there could be large rates for 
gx production, with simple signatures, whereas gg production suffers from large hadronic 
jet backgrounds. 

Experimental investigations are facilitated by firm theoretical understanding of the ex- 
pected sizes of the cross sections for production of the superparticles. In the case of hadron- 
hadron colliders, the large strong coupling strength (as) results in potentially large contri- 
butions to cross sections from terms beyond leading order (LO) in a perturbative quantum 
chromodynamics (QCD) evaluation of the cross section. For accurate theoretical estimates, 
it is necessary to extend the calculations to next-to-leading order (NLO) or beyond. NLO 
contributions generally reduce and stabilize dependence on undetermined parameters such 
as the renormalization and factorization scales. To date, associated production has been 



calculated only in LO , but NLO results exist for hadroproduction of gluinos and squarks 
0, top squarks 0, sleptons |/]||, and gauginos ||. Studies have begun to incorporate these 
NLO results into Monte Carlo simulations 0. 

In this Letter we present the first NLO (in SUSY-QCD) calculation of hadroproduction 
of a g in association with a x, including contributions from virtual loops of colored sparti- 
cles and particles and three-particle final states involving the emission of light real particles. 
We extract the ultraviolet, infrared, and collinear divergences by use of dimensional regu- 
larization and employ standard MS renormalization and mass factorization procedures. In 
the course of computing the virtual contributions, we encountered new divergent four-point 
functions. The contributions from real emission of light particles are treated with a phase 
space slicing method. We provide predictions for inclusive cross sections at Tevatron and 
LHC energies. Our reason to focus on the g + x final state, rather than on the associated 
production q + x, is that at the energy of the Tevatron the LO cross sections for g + x are 3 
to 6 times greater than those for q + x (X. — X121X1) when rrig = rriq = 300 GeV, and 6 to 
15 times greater when m~ g = rriq = 600 GeV. In obtaining the q cross sections, we sum over 
five flavors of squarks and antisquarks. 

In LO of SUSY-QCD, the associated production of a gluino and a gaugino proceeds 
through the subprocess qq — > gx with a t-channel or a w-channel squark exchange. We 
assume that there is no mixing between squarks of different generations and that the squark 
mass eigenstates are aligned with the squark chirality states, equivalent to the assumption 
that the two squarks of a given flavor are degenerate in mass. We ignore the rif = 5 
light quark masses in all of the kinematics and couplings. Under these assumptions, the 
massless incoming quarks and antiquarks have a particular helicity, and thus the Feynman 
diagrams in which a right-handed squark is exchanged cannot interfere with those mediated 
by a left-handed squark. In evaluating the Feynman diagrams involving Majorana and 
explicitly charge-conjugated fermions, we follow the approach of Ref. |10|| . In the case of 
charged gauginos, only the left-handed squarks participate, whereas neutral gauginos receive 
contributions from both left- and right-handed squarks. 



The LO matrix element summed (averaged) over the colors and helicities of the outgoing 
(Incoming) particles has the analytic form || 
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Here, rriq t u is the mass of the squark exchanged in the t- and it-channel; a$ = Osl^R is the 
Yukawa coupling between quarks, squarks, and gluinos. At leading order it is equal to the 
gauge coupling constant as- In Eq. (0) X t> tu,u stands for the Yukawa couplings of quarks, 
squarks, and gauginos. Variables s, t, and u are the usual Mandelstam invariants at the 
partonic level; tg t % = t — m? ^ and Ug ^ = u — fn^^- 

At NLO the cross section receives contributions from virtual loop diagrams and from real 
parton emission diagrams. The virtual contributions arise from the interference of the Born 
amplitudes with the related one-loop amplitudes containing self-energy corrections, vertex 
corrections, and box diagrams. We include the full supersymmetric spectrum of strongly 
interacting particles in the virtual loops, i.e. squarks and gluinos as well as quarks and 
gluons. 

Since the virtual loop contributions are ultraviolet and infrared divergent, we regularize 
the cross section by computing the phase space and matrix elements in n = 4 — 2e dimen- 
sions. We calculate the traces of Dirac matrices using the "naive" 75 scheme in which 75 
anticommutes with all other 7^ matrices. This choice is justified for anomaly-free one-loop 
amplitudes. The 75 matrix enters the calculation through both the quark-squark-gluino and 
quark-squark-gaugino Yukawa couplings. We simplify the integration over the internal loop 
momenta by reducing all tensorial integration kernels to expressions that are only scalar in 
the loop momentum [II]. The resulting one-, two-, three-, and some of the four-point func- 



tions were computed in the context of other physical processes 0. However, we compute 
two previously unknown divergent four-point functions; these new functions arise because 
the final state gluino and gaugino have different masses in general. We evaluate the scalar 
four-point functions by calculating the absorptive parts with Cutkosky cutting rules and the 



real parts with dispersion techniques. 

The ultraviolet (UV) divergences are manifest in the one- and two-point functions as 
poles in 1/e. We remove them by renormalizing the coupling constants in the MS scheme 
at the renormalization scale Q and the masses of the heavy particles (squarks and gluinos) 
in the on-shell scheme. The self-energies for external particles are multiplied by a factor of 
1/2 for proper wave function renormalization. A difficulty arises from the fact that spin-1 
gluons have n — 2 possible polarizations, whereas spin-1/2 gluinos have 2, leading to broken 
supersymmetry in the MS scheme. The simplest procedure to restore supersymmetry is with 
finite shifts in the quark-squark-gluino and quark-squark-gaugino Yukawa couplings [12|,5|. 



In addition to the ultraviolet singularities, the virtual corrections have collinear and 
infrared singularities that show up as 1/e or 1/e 2 poles in the derivatives of the two-point 
function and in the three- and four-point functions. These infrared singularities appear as 
factors times parts of the Born matrix elements. They can be separated into Cp and Nc color 
classes, depending on the color flow and the abelian or non-abelian nature of the correction 
vertices. They are cancelled eventually by corresponding soft and collinear singularities from 
the real three parton final state corrections. 

The real corrections to the production of gluinos and gauginos arise from three parton 
final-state subprocesses in which additional gluons and massless quarks and antiquarks are 
emitted: qq -> gxg, qg -> gxq, and qg -> ~g\q- 

The n-dimensional phase space for 2 — > 3 scattering may be factored into the phase 
space for 2 — > 2 scattering and the phase space for the subsequent decay of one of the two 
final state particles with squared invariant mass s^ = {p\ +P3) 2 —ml into two particles with 
momenta p\ and P3, parametrized in the rest frame of particles 1 and 3 |13| . We follow 



the procedure in Ref. |T3| to reduce all of the angular integrals. The angular integrations 
involving negative powers of t' = {pb—Pzj 2 and u' = (p a ~P3) 2 , where p a and pb are the four- 
momenta of the incoming partons, produce poles in 1/e which correspond to the collinear 
singularities in which particle 3 is collinear with particle a or b. Because these singularities 
have a universal structure, they may be removed from the cross section and absorbed into 



the parton distribution functions according to the usual mass factorization procedure [JT4 . 



In addition to the collinear singularities described above, the corrections involving real 
gluon emission also have infrared (IR) singularities arising when the energy of the emitted 
gluon approaches zero. These singularities appear as poles in S4 in the cross section and 
must also be extracted so that they can be combined with corresponding terms in the virtual 
corrections and shown to cancel. In order to make this cancellation conveniently, we slice 
the gluon emission phase space into hard and soft pieces, 

d 2 a* /-a d 3 a s r* max , d 3 a H 
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where A is an arbitrary cut-off between soft and hard gluon radiation. When the cut- 
off is much smaller than the other invariants, the S4 integration for the soft term becomes 
simple and can be evaluated analytically, leading to explicit logarithms log A/m 2 , log A/m 2 ; 
m = (trig + m^)/2. The hard term is free from infrared and, after mass factorization, also 
collinear singularities and can be evaluated numerically in four dimensions. This procedure 
leads to an implicit logarithmic dependence of the hard term on the cut-off A which cancels 
the explicit logarithmic dependence in the soft term. 

To obtain numerical results for the cross sections, we work within a particular SUGRA 
scheme, but the cross sections depend principally on the masses of the \ and g and are 
otherwise fairly independent of the details of the SUSY breaking. The physical gluino and 
gaugino masses as well as the gaugino mixing matrices are calculated from a default minimal 
SUGRA scenario |l|]. We choose the common scalar and fermion masses at the GUT scale 
to be m = 100 GeV and mu2 = 150 GeV. The trilinear coupling A = 300 GeV, and the 
ratio of the Higgs vacuum expectation values tan/3 = 4. The absolute value of the Higgs 
mass parameter /1 is fixed by electroweak symmetry breaking, and we choose /i > 0. For 
this scenario, we find the neutralino masses m^o to be 55, 104, 283, and 309 GeV with 

Xi_ 4 

m^o < inside a polarization sum. The chargino masses m-± are 102 and 308 GeV and 
therefore almost degenerate with the masses of the m^o and m^o, respectively. 

The total hadronic cross section is obtained from the partonic cross section through the 
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where r = ^§-, and 5 is the square of the hadronic center-of-mass energy (\^S = 2 TeV for 
Run II at the Fermilab pp collider Tevatron and 14 TeV at the CERN LHC pp collider). For 



the NLO predictions, we employ the CTEQ4M parametrization [|J| for the parton densities 
f(x, Q 2 ) in the proton or antiproton and a two- loop approximation for the strong coupling 
constant as with A^ 5 ^ = 202 MeV. In LO we use the LO parton densities CTEQ4L and the 
one-loop approximation for as with A^ 5 ^ = 181 MeV. 

In Fig. p] we present predictions for total hadronic cross sections at the Tevatron (top) 
and the LHC (bottom), as a function of the physical gluino mass. To obtain these results, 
we use the average produced mass as the hard scale Q in Eq.fl3|), Q = (m g + m^)/2. We 
vary the SUGRA parameter m\ii between 100 and 400 GeV and keep the other SUGRA 
parameters fixed to the values listed above. As gluino mass increases over the range shown 
in the figure, the corresponding gaugino mass ranges are 31 to 163 GeV for Xii 62 to 317 
GeV for x 2 an d Xi , 211 to 666 GeV for x°, and 240 to 679 GeV for %° and xf- The chargino 
cross sections are summed over positive and negative charges. 

We observe that the cross sections for x 2 and xf an d those for x\ an d xf are very similar 
in magnitude at the Tevatron, as are their respective masses. One might expect the largest 
cross section for the lightest gaugino Xv However, its coupling is dominantly photino-like 
and smaller than the zino-like coupling of x 2 which therefore has a larger cross section at 
small m g despite its larger mass. The heavier gauginos xt 4 an d Xi are dominantly higgsino- 
like and their cross sections are suppressed by over an order of magnitude with respect to 
those of the lighter gauginos. 

At the Tevatron, the NLO contributions increase the cross sections by 5 to 15% at 
the hard scattering scale Q = (nig + m%)/2, depending on the channel considered and 
the values of the masses. At the LHC, the increases are in the range of 15 to 35%. The 



purely NLO qg incident channel contributes significantly at the LHC, in addition to the qq 
channel, particularly for the lighter gauginos, whereas the qg channel plays little role at the 
Tevatron. In the event sparticles are not observed, the predicted increases translate into 
more restrictive experimental mass limits. 

The enhancements of the cross sections are modest and, as such, underscore the validity 
of perturbative predictions for the processes considered. A further important benefit of the 
NLO computation is the considerable reduction in theoretical uncertainty associated with 
variation of the renormalization and factorization scale Q. For the processes studied here, 
this dependence is typically ±10% at the Tevatron when Q is varied over the interval Q/m 
from 0.5 to 2, compared to ±25% in leading order. At the LHC, the dependences are ±9% 
at NLO and ±12% at LO. 

In this Letter we limit ourselves to total cross sections. Differential distributions in 
the transverse momentum p? and the rapidity r] of the produced sparticles will be pub- 



lished elsewhere [I5fl, along with figures of scaling functions, renormalization/factorization 
scale dependence, K-factors, and several appendices containing a detailed exposition of the 
calculation. 

To summarize, we provide NLO predictions of the cross sections for the associated pro- 
duction of gauginos and gluinos at hadron colliders. If supersymmetry exists at the elec- 
troweak scale, the cross section for this process is expected to be large and observable at the 
Fermilab Tevatron and/or the CERN LHC. It is enhanced by the large color charge of the 
gluino and the (in many SUSY models) small mass of the light gauginos. The cross section 
for x® production is the largest, because of its zino-like coupling. The cross section for xt is 
about equal to that of the x°. The NLO predictions are modestly larger than the LO values 
but considerably more stable theoretically. 

Work in the High Energy Physics Division at Argonne National Laboratory is sup- 
ported by the U.S. Department of Energy, Division of High Energy Physics, under Contract 
W-31-109-ENG-38. The authors are grateful for correspondence with W. Beenakker and 
conversations with S. Mrenna. T. Tait has benefitted from discussions with C.-P. Yuan. 
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FIG. 1. Total hadronic cross sections for the associated production of gluinos and gauginos at 
Run II of the Tevatron (top) and the LHC (bottom). NLO results are shown as solid curves, and 
LO results as dashed curves. We vary the SUGRA scenario as a function of 7774/2 £ [100; 400] GeV 
and provide the cross sections as a function of the physical gluino mass rrig. The chargino cross 
sections are summed over positive and negative chargino charges. 
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